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I.  INTRODUCTION 


The  reaction  zone  in  cw  chemical  lasers  (Fig.  1)  is  generally 
maintained  at  pressures  of  the  order  of  1  to  10  Torr  in  order  to  permit  fast 
mixing  of  the  reactants.  At  these  pressure  levels,  the  spectral  line  shape 
is  inhomogeneously  broadened  (i.  e.  ,  the  radiation  field  interacts  with  only  a 
portion  of  the  excited  molecules).  As  a  result,  hole  burning*  may  occur 
as  the  degree  of  optical  saturation  is  increased. 

2  3 

In  previous  analytical  studies  of  cw  chemical  lasers,  ’  hole -burning 

effects  are  ignored.  In  these  studies,  it  is  assumed  that  lasing  occurs  at 

line  center  of  the  Doppler-broadened  line  shape  and  that  the  latter  line  shape 

is  maintained  in  the  presence  of  lasing.  These  assumptions  provide  for 

reasonable  estimates  of  oscillator  output  power  for  cases  where  resonator 

configurations  permit  a  large  number  of  optical  modes.  However,  for 

4  5 

single-mode  operation  at  pressures  of  the  order  of  1  to  10  Torr  ’  and  for 
multimode  operation  in  relatively  low-pressure  devices,  hole  burning  effects 
need  to  be  considered.  These  effects,  in  fact,  have  been  reported  in  Refs.  4 
and  5,  where  a  "Lamb  dip"  is  observed  as  a  single-mode  cw  chemical  laser 
is  tuned  across  line  center. 

A  comprehensive  theory  for  inhomogeneous  broadening  effects  in  a 
steady-state  laser  oscillator  has  been  developed  by  Lamb.  ^  Lamb's  theory 


is  directly  applicable  to  atomic  lasing  systems  wherein  the  deactivation  of 
excited  particles  (by  spontaneous  emission)  is  fast  compared  to  the  particle 


collision  rate  (i.  e.  ,  the  effect  of  collisions  on  the  excited  particle  velocity 
distribution  function  is  neglected).  Kan  and  Wolga  extended  Lamb's  theory 

to  CO  molecular  lasers  wherein  the  deactivation  rate  of  excited  CO  is 

6  2 

slow  compared  to  the  molecular  collision  rate.  The  latter  collisions  tend 
to  fill  the  "hole"  induced  by  inhomogeneous  broadening  and  need  to  be 
considered. 

G  7 

In  the  previous  theories  ’  for  inhomogeneous  broadening  effects, 
flow  properties  do  not  vary  temporally  or  spatially.  The  extension  to  cw 
chemical  lasers  is  not  straightforward  since  flow  conditions  in  the  latter 
are  functions  of  streamwise  distance.  Hence,  the  present  study  was 
undertaken  to  evaluate  the  effects  of  inhomogeneous  broadening  on  cw 
chemical  laser  performance.  A  simple  two-level  vibrational  model 
similar  to  that  used  in  Ref.  2  is  used.  A  Fabry-Perot  (F-P)  resonator  is 
assumed,  and  the  variation  of  flow  conditions  with  streamwise  distance  is 
considered.  Discrete  rotational  energy  states,  however,  are  not  considered. 


II.  THEORY 


The  flow  is  considered  to  have  a  uniform  velocity  u  in  the  +x  direction. 
A  Fabry-Perot  (F-P)  resonator  is  aligned  with  its  optical  axis  normal  to  the 
flow  (Fig.  1).  The  total  number  of  particles,  per  unit  volume,  in  the  lower 
and  upper  vibrational  levels  are  denoted  n^  and  n^>  respectively.  General 
aspects  of  the  interaction  between  the  radiation  field  and  the  lasing  medium 
are  discussed.  Effects  of  inhomogeneous  broadening  on  the  performance  of 
a  cw  chemical  laser  are  then  deduced.  Steady-state  lasing  is  assumed. 


A. 


GENERAL  CONSIDERATIONS 


In  the  following  sections,  inhomogeneous  broadening  and  the  longitudi¬ 
nal  mode  structure  in  a  F-P  resonator  are  discussed. 

1.  STATIONARY  PARTICLES 

Let  Vq  denote  the  resonant  (line  center)  frequency  for  absorption  or 
stimulated  emission  of  radiation  by  particles  essentially  at  rest.  (More 
specifically,  collisions  are  permitted  between  particles,  but  terms  of  order 
v  /c  compared  to  one  are  neglected.  )  For  these  particles,  the  gain  at 

y 

frequency  v  can  be  expressed 


g(v,vQ)  =  o(v,vQ)  (n^  -  n^ ) 


where  o(v,Vq)  is  the  cross  section  for  stimulated  emission  at  v.  Equation  (1) 


can  be  written  in  the  form 


8 


g(v,vQ)  =  [^(v  -  v0)]  (n2  -  n^ 


(2a) 
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where  =  ct  (vq,  Vq)  is  the  cross  section  at  Vq  and  £?(v  -  Vq)  is  the 
Lorentzian  (homogeneous)  line  shape 


U’(v  -  vQ)  = 


CT(v,  V/J 


Here,  Av^  is  the  characteristic  width  (FWHM)  of  the  homogeneous  line  shape 
and  is  a  function  of  pressure  (Appendix  A),  The  broadening  is  due  to  particle 
collisions.  The  index  of  refraction  at  v  can  be  expressed 


TKv,v0)  -  1  =J^  g(v’V 


which  follows  from  the  Kramers- Kronig  relations. 

The  integral  of  the  homogeneous  line  shape  is  a  constant.  Thus, 


a(v,  ,Q)  dv  =  ^  aQ  Avh 


where  Av^  is  independent  of  pressure  level,  but  Av^  and  a  are 
proportional  to  pressure. 

2.  DOPPLER  EFFECT 

The  effect  of  random  motion  of  the  particles  is  considered.  1  +(v) 

denotes  radiation  with  frequency  v  traveling  in  the  +y  direction,  and  v  the 

random  particle  velocity  in  +y  direction  (Fig.  1).  Because  of  the  Doppler 

effect,  the  frequency  v>  which  will  be  resonant  with  particles  with  a  velocity 

v  is  v  =  v  /f  1  -  (v  / c ) ]  which,  for  the  realistic  assumption  v  /c«  1,  becomes 
y  o  y  y 
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V 


(3a) 


1  + 


Conversely,  the  particle  velocity  v  ,  which  will  result  in  resonance  with 
I  +(v),  is 


V 

- 

C  “ 


fe)- 


(3b) 


In  the  absence  of  radiation,  the  particles  in  the  upper  and  lower  vibrational 
levels  can  be  assumed  to  have  a  Maxwellian  velocity  distribution.  The  num¬ 
ber  of  particles  resonant  with  frequencies  in  the  range  v  to  v  +  dv  is  then 
n(v)  dv,  where^ 


n(v)  _  1 

n  Av 


14  In  Z\i,Z 

\T~J  exP 


2 1 


(4  In  2) 


Avd 


(4) 


Here,  Avd  is  the  characteristic  (Doppler)  width  (FWHM)  of  the  Maxwellian 
distribution  [Eq.  (A-  1 )].  For  later  use,  the  notations 


are  introduced.  Note  that 


p(v)  dv 


p(v)  dv  =  1 


(5c) 


Upon  reflection  from  a  lossless  mirror,  I  +(v)  will  be  converted  to  7~(v), 
which,  as  a  result  of  the  Doppler  effect,  interacts  with  particles  with 
velocity  v^/ c  =  -[(v)/vq)  -  1].  Thus,  I  ~  (v)  interact  with  the  velocity  groups 
vy/c  =  ±[(v/vQ)  -  1],  respectively.  Note  also  that  I  ‘(v)  andT  +  (2v0  -  v) 
interact  with  the  same  particles.  This  equivalence  is  used  in  the  subsequent 
treatment  of  an  F-P  resonator.  In  particular,  particle  groups  are  identified 
by  the  resonant  frequency  associated  with  the  wave  I  +(v).  Also,  integration 
over  all  population  groups  is  accomplished  by  integration  with  respect  to  v 
in  the  interval  <  v  <  05  rather  than  by  integration  with  respect  to  v  in 

y 

the  interval  -°°  <  v  <  «>. 

y 

3.  INHOMOGENEOUS  BROADENING 

Under  lasing  conditions,  the  particle  distribution  function  is  perturbed 
from  the  Maxwellian,  and  more  general  distributions  need  to  be  considered. 
Thus,  from  Eq.  (2),  the  gain  and  refractive  index  at  v'  resulting  from 
particles  that  are  resonant  with  radiation  in  the  range  v  to  v  +  dv  can  be 
expressed  as  (Fig,  2) 
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tractive  Index  at  v'  as  Result  of 
l  Frequency  Interval  v  to  v  +  dv 
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The  net  gain  and  refractive  index  at  v  '  is  then  found  by  integration  over  all 
parti  le  groups.  The  result  is 


/® 

X’(s> '  -  v)  [n^  (v)  -  n^  (v)]  dv 


T)(v')  - 1  =~z^f  ^r-Av' " v)  [nz(v) "  nl  (v)]  dN 


Equation  (7)  defines  the  gain  and  refractive  index  of  an  inhomogeneously 

broadened  medium.  The  evaluation  of  Eq.  (7)  for  a  chemical  laser  medium 

— ±  / 

is  the  subject  of  the  present  study.  Amplification  of  I  (v  )  is  found  from 


dl  ~  ( v 


=  ±  g(v  )  I  -(v') 


For  a  nonlasing  medium,  Eq.  (4)  is  applicable,  and  Eq.  ,  ’)  be  omes 


(W2)  [a 


/•QD 

- - -  /  p(v)  !£( v 1  -  v) 

0Po  Avh  (n2  -  nj)]  -  *4  ^  Avh 


(4ir/V 
(tt/2)  [aQp 


v  r^Kv r>  -  n  _  4  /*“  v/ 

PoAvh  (n2  -  nl)]  *  L  * 


~  P(v)  S?(v'  -  V) 

Avh  Avh 


The  right-hand  sides  (RHS)  of  Eqs.  (9a)  and  (9b)  equal  the  real  and  imaginary 
part  of  a  complex  error  function.^’  ^  For  (p^  Av^)^  «  1,  fV,°se  quantities 
can  be  deduced  from  Eqs.  (B-ld)  and  (B-le)  with  0^  =  1  therein. 
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4. 


RESONATOR  MODES 


An  F-P  resonator  can  support  longitudinal  modes  centered  at  frequen¬ 
cies  y  =  Nc/2E,  where  N  is  an  integer  and  L  is  the  separation  between 
mirrors.  The  longitudinal  mode  separation  is  then  Avc  =  c/2L.  A  typical 
longitudinal  mode  distribution  is  illustrated  by  the  solid  lines  in  Fig.  3.  The 
frequency  width  of  each  mode  can  be  assumed  to  be  small  compared  to 

Av 

A  convenient  notation  for  dealing  with  the  interaction  between  the 
resonator  modes  and  the  gain  medium  is  now  introduced.  Without  loss  of 
generality  it  can  be  assumed  that  one  longitudinal  mode,  centered  at  Vj,  is 
located  in  the  interval 


v0  < 


vi< 


(10) 


The  central  frequency  of  the  longitudinal  modes  can  be  expressed  in  the  form 


j  =  1,  3,  5,  ... 

j  —  -2,-4,-b,  ... 


(11a) 

(lib) 


These  frequencies  are  denoted  by  solid  lines  in  Fig.  3.  The  radiation 

-  +  t 

I  (vj)  interacts  with  particles  characterized  by  the  frequency  v  =  Vj  in 

Fig.  3.  The  interaction  results  in  the  "holes"  indicated  therein.  As  pre¬ 


viously  noted,  the  interaction  of  the  reflected  radiation  I  (v  )  with  the 
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Fig.  3.  Longitudinal  Mode  Structure  in  F-P  Resonator. 

Solid  lines  denote  longitudinal  modes.  Dashed 
lines  denote  additional  r  odes  for  equivalent 
resonator  with  radiation  in  +y  direction  only. 
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lasing  medium  is  equivalent  to  the  interaction  of  I  +  (2vq  -  v.)  with  this 
medium.  An  equivalent  I  (vT  radiation  field  with  frequencies  defined  by 

vj  =  Zv0  '  V1  +  J  (-~r)  j  =2,4,6,...  (12a) 

=  2V0  -  +  (j  +  1)  (— ^)  j  =  -1,-3, -5,  ...  (12b) 

is  considered.  The  latter  are  indicated  by  the  dashed  lines  in  Fig.  3.  These 
lines  are  the  mirroi  image  of  the  solid  lines,  about  Vq,  and  also  have  a 
separation  Av^.  The  results  of  the  present  approach  are  identical  with  those 
obtained  from  an  F-P  resonator  with  both  I  +  and  1  "  radiation. 

A  further  simplification  can  be  introduced.  Because  of  the  symmetry 
about  Vq,  only  lasing  frequencies  in  the  interval  >  Vq  need  to  be 
considered.  These  are  given  by  Eqs.  (11a)  and  (12a).  In  particular, 
j  =  1,2,3,  ...  is  considered,  where  is  the  value  of  j  for  the  last  (final) 

lasing  transition  in  the  interval  Vj  >  Vq.  The  total  number  of  lasing 
transitions  is  2jf.  The  sum  of  a  quantity  (  )^  over  all  j  is  then 

j=±jf  j£ 

’j  *£  "j  ■  2E '  ’j  ,ni 

J  j=±i  j=i 


When  lasing  transitions  are  separated  by  a  frequency  difference  of  order 

Av,  ,  they  compete  for  the  same  particles.  Later,  Av,  «  Av  is  assumed 
n  he 


so  that,  at  most,  only  two  radiation  fields  can  compete.  Two  cases  can  be 


distinguished.  When  0  <  v£  -  Vq  <  (A^/ 4),  the  separation  between  adjacent 
modes  Av  is  (Fig.  3) 


Avs  =  2  (vt  -  vQ) 


(14a) 


If  Av  /av,  =0(1),  competition  will  occur  between  v  ,  and  v,  ,  v,  and  v0, 
s  n  - 1  1  l  3 

v4  and  v5,  .  .  .  .  Here,  jf  is  odd.  When  (Avc/4)  <  -  vQ  <  (A\>J2.),  the 

separation  between  adjacent  modes  is  (Fig.  3) 


Av 

s 


V1 


(14b) 


If  Avg/Av^  =  0(1),  compeunun  occurs  between  v£  and  v^  and  v^,  .  .  .  . 

Here,  j£  is  even.  When  Avg»Av£i,  no  mode  competition  occurs. 

The  open  interval  is  used  in  Eq.  (10;  to  avoid  the  degenerate  case 

where  two  lasing  transitions  overlap  in  Fig.  3  (i.  e.  ,  Av  =  0).  This  case 

s 

is  treated  herein  by  consideration  of  the  limit  Av  — •  0. 

s 

5.  RESONATOR  BOUNDARY  CONDITION 

Diffraction  effects  are  neglected  and  it  is  assumed  that  each  mirror 
of  the  F-P  resonator  has  the  same  reflectivity  R^.  Later,  it  is  assumed 
that  there  are  ngc  semichannels,  each  of  width  w,  and  that  the  lateral  width 
of  the  gain  region  in  each  semichannel  is  y£(x)  (Fig.  4).  Under  steady-state 
lasing  conditions,  the  net  gain  per  pass  equals  the  net  loss  per  pass  at  each 
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streamwise  station  where  lasing  occurs.  For  these  conditions,  the  local 


gain  is  defined  by 


-In  R 

_ m_  _ 

"j  ~  yf(x)  =  gc 


whe re  g.  =  g(v. ). 

J  J 

B.  GENERAL  SOLUTION 

The  equations  that  define  the  effects  of  gain  saturation  on  the  perfor¬ 
mance  of  a  Doppler-broadened  cw  chemical  laser  are  deduced.  A  premixed 
flow  is  first  discussed.  Effects  of  diffusion  then  are  introduced. 

1.  PREMIXED  LASER 

At  sufficiently  low  pressures,  the  rate  of  diffusion  can  be  considered 
fast,  relative  to  the  rate  of  the  chemical  reactions  in  a  cw  chemical  laser. 
For  these  conditions,  the  reactants  can  be  considered  to  be  premixed  and 
to  start  reacting  at  x  =  0.  Fluid  properties  are  a  function  only  of  streamwise 
distance  x.  The  rate  of  change  of  n2(v)  and  nj(v)  with  x  can  then  be 
expressed 


u— —  =  p(v)  u  -  kcd  n2(v)  +  k__  [p(v)  n,  -  n,  (v)] 


n2  "  nl^ 


1'  V"  - 

—  E  V'Y 


dn^(v)  _ 

U  dx  =  °  +  kcd  n2  (v)  +  kcr  [p(v)  rtj  -  nj  (v)] 


n,  (v)  -  n.  (v) 


A  \  V/  ^ ^  — 

— £  W'" 


(16b) 


where  I.  si  +(v.)  and  n  =  n,  +  n?.  The  terms  on  the  right-hand  side  of 

Eq,  (16a)  have  the  following  interpretation.  In  the  first  term,  the  chemical 

pumping  reaction  (which  is  considered  known)*"  is  assumed  to  create  only 

n^  particles,  which  initially  have  a  Maxwellian  velocity  distribution.  The 

second  term  represents  the  loss  resulting  from  collisional  deactivation.  The 

rate  c  efficient  k^  is  assumed  to  be  independent  of  v.  The  third  term 

represents  the  creation  of  n^(v)  as  a  result  of  cross  relaxation.  The  latter 

is  assumed  equal  to  the  product  of  a  cross  relaxation  rate  coefficient  k 

with  the  difference  between  the  equilibrium  Maxwellian  value  and  the  local 

value  of  n^(v).  This  model  is  similar  to  that  used  in  Ref.  11,  in  which  a 

steady- state  four-level  amplifier  was  treated.  The  last  term  represents 

the  loss  of  n  (v)  as  a  result  of  stimulated  emission  centered  v-.  This 
^  J 

term  is  deduced  from  Eq.  (8)  and  is  nonzero  only  when  v  -  Vj  iO(A^). 

When  n  is  in  moles  per  unit  volume,  the  quantity  e  equals  the  energy  per 
mole  of  photons.  The  terms  on  the  right-hand  side  of  Eq.  (16b)  can  be 
deduced  from  the  preceding  discussion. 

2.  DIFFboION  EFFECT 

A  simplified  diffusion  model  is  adopted  (Fig.  4).  It  is  assumed^’  ^ 
that  the  chemical  reactants  are  premixed,  but  do  not  start  to  react  until 
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' 


Fig.  4.  Simplified  Model  of  CW  Chemic  1  Laser. 
A  single  semichannel  is  shown. 


they  enter  the  reaction  zone  bounded  by  yf  =  yf(x).  Laminar  and  turbulent 
mixing  are  modeled  by 


(17) 


where  m  -  1/2  and  1  for  laminar  and  turbulent  flow,  respectively.  Here, 
w  is  the  semiwidth  of  a  single  oxidizer  (F)  nozzle,  and  xD  is  the  distance 
at  which  the  reaction  zone  reaches  the  nozzle  centerline.  The  latter  is  a 
measure  of  diffusion  rate  and  is  assumed  to  occur  downstream  of  the  lasing 

region. 

When  mean  rates  are  used,  flu:u  properties  within  the  reaction  zone 

may  be  assumed  to  vary  only  with  x.  For  species  n.,  the  variation  can  be 

,12,  13 
expressed 


u 

yf 


d[(n.  -  (n.)J  yf 
dx 


w. 

l 


(18) 


where  (n.)  is  the  number  density  upstream  of  the  flame  sheet,  and  w  is 

'  i  oo  1 

the  volumetric  production  rate  of  m  resulting  from  chemical,  collisional,  and 
radiative  processes.  Note  that  [n^v)]^  =  0;  Eq.  (16a)  becomes 


u_  dtn2(v)yf]  up(v) 

y£  dx  yf 


d(nT  yf) 
dx 


cd 


n2(v) 


+  kcr  [pM n2 


n2(v)] 


n2  "  nl 

€ 


S  a(vj’ v)  h 

i 


(19) 
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A  similar  equation  is  deduced  for  n^(v).  The  following  normalized  unit 
order  variables  are  introduced 


n2  yf 
N  =  - 

2  n  w 


N 


n2(v)_yf 

2v  n  w  p” 
r  *o 


n„  v- 

T  J'f 

N  =  — - — — 
T  n  w 


N 


yf 

1  v  n  w  'p 
r  c 


k  ,  x 
cd 

u 


’D 


k  J  XT^ 

cd  D 
u 


<£. 

J 


v) 


I.  = 
J 


I  . 


e  k 


cd 


(20a) 

(20b) 

(20c) 

(20d) 

(20e) 


G.  =  - 
J  o 


hlL 


-In  R 


m 


n  n  Av,  p  w 
Or  h  ro 


Jh  Fo 


sc 


(20f) 


Here,  n  is  a  characteristic  number  density,  which  for  a  chemical  laser  is 
r 

taken  to  be  equal  to  the  F  atom  concentration  upstream  of  the  flame  sheet 

(n^  =  [F]^).  The  local  lasing  intensity  F  is  normalized  herein  by  e  ^ccj/ao* 

Thus,  the  quantity  L  is  the  ratio  of  the  stimulated  emission  rate  to  the  rate 

of  collisional  deactivation  of  particles  in  the  frequency  interval  about 

v..  An  alternate  normalization,  not  used  herein,  is  1=1/1  where 
J  s 

I  =  [(1  +  R ) / 2]  e  ^c(j/aQ  the  appropriate  measure  of  saturation  intensity 
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in  the  present  model  [e.  g.  ,  Eq.  (28)].  Thus  21/(1  +  R)  is  a  measure  of 
degree  of  saturation  in  the  present  model.  The  normalized  gain  at  v.  is 
obtained  from 


- 


V(N2v-N1v> 


dy 

Avh 


(21a) 


In  the  absence  of  lasing,  N-,^  -  =  p(v)  (N^  -  Nj),  and 


G. 

J 


(N,  - 


N, 


p(v)y'(v 


V. 

J 


dv 

A\ 


(21b) 


When  Av^/Av^  «  1,  Eq.  (21b)  becomes 


G. 

J 


VN2 


Nl> 


(21c) 


where  p.  =  p(v.).  Threshold  is  reached  when  G.  =  G  . 

J  J  J  c 

Substitution  of  Eq.  (20)  into  Eq.  (19)  and  the  equivalent  expression  for 
rij(v)  yields 


d  N?  dN 

ST  "  P,V)  -dT  '  N2v+  R  [PM  N2  '  N2v] 


-  <"2.-NK’£Vi 


(22a) 


(22b) 


d  N1 

“dT  =  N2v  tR[plv>  N1  -  N1V]  *  (N2v  -  N1v>E  Vj 

j 

where  d  NT/d£  is  known.  Integrals  of  Eq.  (22)  form  the  basis  of  the  present 
study.  Addition  of  Eqs.  (22a)  and  (22b)  and  integration  with  respect  to  Q 
and  v  yields,  respectively, 

N2v+  Nlv  =  P(v)  NT  (23a) 

and 

N2  +  Nt  =  Nt  (23b) 

The  fact  that  Eq.  (23a)  is  independent  of  R  is  somewhat  surprising  and  indi¬ 
cates  that  the  increase  in  N^,  as  a  result  of  cross  relaxation,  is  offset  by 
an  equivalent  loss  in  N^.  Equation  (23b)  is  consistent  with  the  definition  of 
N-p.  Subtraction  of  Eq.  (22b)  from  Eq.  (22a)  and  integration  with  respect 
to  v  yields,  respectively, 

d  (N  -  N.  )  [dN 

- - -  +  <i  +  R)  (n2v-  n1v)  =  p<v)|-jU  .  Nt  +  r,^  .  Nj) 

- 2  <n2v  -  Niv’z: ^ 

j 
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and 


d  (N  -  N  )  dN 

— 3c — *  (N2  -  Ni>  =-ar  -  nt  -  2<p0AV  Gc‘  (24bl 


The  last  term  on  the  right-hand  side  of  Eq.  (24b)  follows  because  G  =  G 

J  c 

when  I  *  0.  The  zero  power  solution  I  =  0  of  Eq.  (24b)  is 


which,  together  with  the  relation  =  N^, ,  defines  and  N^.  In 

order  to  evaluate  Eq.  (24)  for  I  *  0,  note  that  for  Ij  *  0  (i.  e.  , 

Gj  -  Gc>’ 


Equation  (26a)  is  satisfied  by  d  (N?  -  N,  )/d£  =  0,  which  is  a  sufficient,  but 

^  V  1  v 

not  a  necessary,  condition,  suggesting  that  for  I.  *  0 


1  (26b) 


-28- 


be  assumed  in  Eq.  (24a).  The  latter  becomes,  with  j  used  as  the  summa 
tion  index. 


N2v  '  Nlv 


P(v)[(dNx/dC)  -  Nt  +  R(N2  -  Nt)] 

1  +  R  +  2  ^2  V(v  -  v./)  I.  / 

.  /  J  J 

J 


Multiplication  by !f(\>  -  v.)  dv/Avh  and  integration  yields,  for  I.  *  0, 

J  J 


p(v)  f/'(v  -  \>-)  dv/Av, 


Pj  (dNT/dC)  -  Nt  +  R(N2  -  Nt)  Pj  1  +  [2/(1  +  R)]  -  v./)  E 


Equation  (28)  provides  implicit  relations  between  L  and  N2  -  N^f  which 
can  be  used  to  integrate  Eq.  (24b).  This  integration  provides  N2  -  and 
I.  as  functions  of  Q,  and  all  other  laser  properties  of  interest  can  be  readily 
deduced.  For  example,  the  power  emitted  in  the  interval  0  S  x^  <  x  by  a 
laser  of  unit  height  is  found  from 


(x)  =  nsc  f  ^fgc 


The  validity  of  Eq.  (26b)  for  cw  chemical  lasers  can  be  further  supported 
by  noting  that  all  quantities  in  this  equation  are  of  order  one  except  for  R 
which  can  be  considered  large  (Appendix  A).  The  validity  of  Eq.  (26b)  for 
moderate  values  of  R,  however,  requires  further  study. 


1 


In  nondimensional  variables,  Eq.  (29a)  becomes 


P(0 


(29b) 


The  streamwise  station  at  which  lasing  ends  (I  =  0)  is  denoted  £e  and  the 
cor  responding  laser  net  output  power  is  denoted  P  .  Similarly,  the  refrac¬ 
tive  index  in  the  lasing  region  can  be  obtained  from 


1  +  R  2tt  [T1(V  -  1]yf/w 
X 


J 


o_  n  pn  Av, 
0  r  0  h 


dN 


T 


dC 


-  nt  +  r(n2  - 


Nil] 


ifi 

Pi  L  ' 


j  p(v)  y(v.  -  v)  dv/Avh 

^Vh  i  +[2/(1  +  R)]]T  -  v)  I. 


(29c) 


For  cases  with  no  mode  competition,  Avg  »  Av^,  the  summation  in  Eqs.  (28) 
and  (29c)  contains  only  a  single  term  and  the  integrals  can  be  evaluated  in 
closed  form  (Appendix  B). 

The  solution,  of  the  system  of  equations  defined  by  Eqs.  (24b)  and  (28) 
generally  requires  numerical  integration.  Closed  form  solutions  can  be 
obtained  in  limiting  cases,  however,  and  these  are  discussed  in  the  next 
section. 

C.  LIMITING  CASES 

The  assumptions  Av^  <<  Av£  and  Av^  «  Av^  are  introduced  in  order 

to  simplify  the  integration  of  Eq.  (28).  The  limiting  cases  [2L/(1  +  R)]^  «  1 

and  Av  «  Av,  are  then  discussed, 
s  n 
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When  Av^«  Av^,  at  most,  two  radiation  fields  compete  for  the  same 
particles.  These  fields  are  denoted  j  '  =  j  and  j '  =  j  +  1  in  Eq.  (28). 

From  the  assumption  Av^«  Av^,  it  follows  that  p(v)  =  p j ,  and,  when  two 
fields  compete,  T  -  L  +  j.  The  right-hand  side  (RHS)  of  Eq.  (28)  then  becomes 


RHS 


lf(\>  -  v.)  dv/A 
1  +  [2/(1  +  R ) J  I.  L'/'(v 


v.)  +S/U  -  vj+1)J 


J 


(30) 


where,  from  Eq.  (14),  v- ,  ,  -  v.  =  Av  .  Two  limiting  cases  are  treated, 

j  +1  j  s  ° 

2 

i.  e.  ,  [21.  /  ( 1  +  R)]  «  1  and  Av  «  Av,  . 

S  il 

1.  CASE  [ 2T / ( 1  +  R)]2  «  1 

The  present  case  corresponds  to  weak  saturation.  With  the  expansion 
(1+e)  *=l-e+0  (e2)  and  the  use  of  Appendix  B, 


RHS 


=  1  )l  .  -h 

2  I 


1  +  R 


1 1  '(t1) 


•  (All 


(31) 


Substitution  into  Eq.  (28),  taking  the  reciprocal  of  both  sides,  and  solving 


for  L  yields, 


21. 

_ 1_ 

1  +  R 


7T 

2 


[(dNx/d£)  -  Nt  +  R(N2  -  Nt)]  p. 

(1  +  R)  G  ‘  1 

c 


(32) 


The  total  intensity  is,  from  Eq.  (13), 
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21 

1  +■  R 


w  r(dNT/d£)  -  Nt  +  R(N2  -  Nj)]^ 

2  (i  +  R)  Cc  Pj 

j  =  l 


h 

(33) 


The  present  solution  is  self-consistent  when  the  square  of  the  right-hand 
side  of  Eq.  (32)  is  small  relative  to  one. 

A  closed-form  solution  of  Eq.  (24b)  can  be  obtained  in  the  present 
limit.  With  the  constants 


A  =  1  + 


(34a) 


B  =  1 


(34b) 


C  = 


4jf  (1  +  R)  Gc  pQ  Avh 


1  +  </■ 


introduced,  Eqs.  (33)  and  (24b)  become 


(34c) 


2  (Pn  Av,  )  G  I 


0  c  A  “  R  \  dC 


I^X  .  N  \ 

R  \  dC  ntJ 


nJ  +  (A  -  1)  (N2  -  Nt)  -  C  (35) 


d(N2  -  Nj) 

dC 


+  A (N2  -  N 


i>  ■ B  (^r  ‘  nt)  * c 


(36) 
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Equation  (36)  applies  to  the  nonlasing  as  well  as  the  lasing  regime.  Lasing 
is  initiated  on  a  given  transition  Vj  at  the  station,  where  G.  first  becomes 
equal  to  Gc  [Eq.  (21c)].  At  small  Q,  -  Nj  is  small  and  threshold  is  not 
reached.  Here,  A  =  B  =  1,  C  =  0,  and  Eq.  (36)  agrees  with  Eq.  (24b),  with 
1  =  0.  Lasing  is  first  achieved,  by  definition,  on  the  Vj  transition.  The 
quantities  A,  B,  and  C  are  then  re-evaluated  with  =  1  and  remain  con¬ 
stant  until  the  second  transition  is  initiated.  Thus,  A,  B,  and  C  are  piece- 
wise  uniform  in  intervals  £.  <  G  <  Ci+1»  which  are  defined  by  the  initiation 
or  termination  of  lasing  transitions.  The  piecewise  integral  of  Eq.  (36)  for 
the  interval  £.  to  +  is 


I  ,N2  ‘  N 


AC.  fC  AC  dN 

.)  e  0  =  (A  +  1)  /  e  dG 

1  J  C.  J  C.  0  0 


(b  -  nt) 


Acor 


The  net  power  emitted  in  this  interval  is  from  Eqs.  (29)  and  (35) 


A  -  1 


;0.f  =(NT-7T?rco)C  +jf[  R(N2  -  N,)  -  NT]  d£0 

Jc,  '  C,  Ai 
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where  is  obtained  from  Eq.  (38).  Equations  (23),  (35),  (38),  and 

(39)  define  the  performance  of  a  chemical  laser,  provided  is  known. 

2.  CASE  Av  «  Av, 
s  h 

In  this  case,  the  competing  lasing  modes  have  nearly  the  same 

frequency.  Thus,  1\\>  -  v.  . )  =j/,(v  -  v.)  [1  +  0  (Av  /Av,  )]. 

J  '  i  J  s  n 

Integration  of  Eq.  (28)  and  the  use  of  Appendix  B  yields 


2  1. 

_ 1_ 

1  +  R 


1  ( 


[(dNT/dC)  -  Nt  +  F.' 
2  (1  +  R)  C 


2iZi 


J 


(40) 


and 


2  I 

1  +  R 


(dNx/dC)  -  Nt  +  R(N2  - 


N. 


(1  +  R)  G 


jf 

E 

j=l 


2  . 
P;  '  h 


(41) 


In  the  limit,  [ I. / ( 1  +  R)]^  «  1,  Av  — 0.  Equations  (32)  and  (40)  are  in 
J  S 

agreement,  as  expected.  Substitution  of  Eq.  (41)  into  Eq.  (24b)  yields  a 
first-order  nonlinear  equation  that  requires  numerical  integration. 

D.  PUMPING  RATE 

An  expression  for  is  needed.  A  useful  limiting  case  is  to  assume 
that  the  chemical  reaction  that  creates  excited  species  (i.  e.  ,  the  pumping 
reaction)  is  fast  relative  to  the  diffusion  rate.  In  the  present  model,2’  *2 
this  reaction  is  equivalent  to  the  assumption  that  (i.  e.  ,  fp]  )  is  con- 
verted  to  n^  (i.  e.  ,  [HF  ])  at  the  instant  that  the  n^  particles  enter  the 
reaction  zone  bounded  by  y^.  It  follows  that 
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N 


T 


(42) 


where  m  =  1/2  and  1  for  laminar  and  turbulent  flows,  respectively. 

E.  LAMINAR  DIFFUSION 

Laminar  diffusion  is  considered  and  analytic  solutions  deduced  for  the 
zero-power  case  and  for  [2L/(1  +  R)]  «  1. 

1.  ZERO  POWER 

Substitution  of  Eq.  (42),  with  m  -  1/2,  into  Eq.  (25)  yields 


-  1/2 


(n2  -  Ni> 


2D(C1/Z)  - 


r  1  /2 

9 


(43) 


where  D(  )  is  the  Dawson  integral.  Since  gain  saturation  effects  do  not 
enter  into  the  zero-power  solution,  the  present  results  are  identical  to  those 
of  Ref.  2.  The  population  inversion  has  a  maximum  £^^2  (N^  -  Nj)  =  0.3528 
at  £  =  0.  3051.  The  inversion  is  zero  at  £  =  1.  1301.  Further  results  are 
given  in  Ref.  2.  The  station  at  which  lasing  is  initiated  can  be  determined 
from  Equations  (43)  and  (21).  Thus,  for  j  =  1,  lasing  is  initiated  at  the 
station  £^,  which  is  defined  by 


2 

IT 


’D 


1/2  (N2  -  Nj).  =  2D(£.1/2)  -  £.1/2 


(44) 


from  which,  it  follows  that  0  <  £^  <  0.  3051.  Values  of  that  permit 
lasing  are  defined  by  2 *  ^2  /irp  ^  <  0.  3528. 
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2. 


CASE  [  2  L  /  ( 1  +  R)]2  «  1 

Substitution  of  Eq.  (42)  into  Eq.  (38)  yields 


A 
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Cd2  [(N2  -  N,> 


A  +  1 
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(45a) 
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dC 


0=(' 


A  +  1  Ll/2 
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C  1/2 
B 


;-A(C  -  C.)  _  j 


D[<Ah>1/2]  -  c‘/z 


C  -1/2_A  -1/2.  _N  ) 

+  B  B  '  2  l’l 


(45b) 


Substitution  of  these  expressions  into  Eqs.  (35)  and  (39)  yields  explicit 
expressions  for  the  local  lasing  intensity  and  net  output  power.  It  is  of 
interest  to  evaluate  the  lasing  intensity  at  the  station  where  lasing  is 
initiated.  Substitution  of  Eq.  (35)  [note  Eq.  (34)]  yields  (with  p^  =  pt  and 

C  =Ci> 


(46) 
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_  V  Ji 


2 

For  Q  near  0.  3051,  I  is  small  and  the  assumption  [1/(1  +  R)]  «  1  is 

satisfied  for  all  values  of  R.  For  small  the  value  of  R  must  be  large 
in  order  for  the  solution  to  be  self-consistent.  The  present  solution  is  now 
examined  in  the  limit  Q  -  small  and  in  the  limit  R  —  °°. 

Consider  Q  -  small  and  R  *  <».  The  net  power  emitted  up  to  station 
C  is  found  by  expanding  I  in  a  Taylor  series  about  and  substituting  into 
Eq.  29b.  The  result  is 


p(C)  =  PQ  Avh  Gc  I. 


(C 


Ci> 


!i 


+  o 


7  dl  \ 

\dq 


(C  -  C :  ) 


(47) 


Thus,  the  output  power  is  proportional  to  p^  and  is  independent  of  R. 
In  this  case,  particle  residence  time  in  the  lasing  region  is  too  short  for 
cross  relaxation  to  contribute  to  the  laser  output. 

In  the  limit  R  — °°,  the  solution  becomes,  for  =  1, 


N. 


N1  = (N2  - 


Ni>i  + 


O 


(a.) 


(48a) 


2  (pn  AvJ  CV2  C  I  = 


2"pl  (P0  AV 


1  + 


—  r 1  -  C1  (N  -  N  ) 

TJI  fei  u  2  l'i 


for  C  =  C: 


(48b) 


(Continued) 


In  this  limit  the  cross  relaxation  is  sufficiently  fast  that  the  lasing  medium 

acts  like  a  homogeneous  medium.  There  is  no  hole  burning,  and  the  results 

are  the  same  as  those  deduced  in  Ref.  2.  Note  that  I  is  discontinuous  at 

C  =  Cj.  a  result  of  the  assumption  R  — -  and  the  corresponding  neglect  of  the 

term  of  order  e  in  Eq.  (48c). 

3.  CASE  Av  «  Av, 
s  n 

Laminar  flow  solutions  are  obtained  by  the  substitution  of 

1/2 

N-p  =  into  Eqs.  (41),  (24b)  and  (29).  Numerical  integration  is 

required. 
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HI.  RESULTS  AND  DISCUSSION 


The  parameters  deduced  herein  are  similar  to  those  in  Ref.  2,  except 

for  the  additional  parameters  pQ  Av^  =  0.  9394  and  R.  Typical 

values  for  cw  chemical  lasers  are  discussed  in  Appendix  A.  In  particular, 

R  =  O  (100)  and  Av^/Av^  =  O[0.0i  X  p(Torr)],  where  p  is  the  net  static 

pressure  in  the  lasing  region  and  generally  is  in  the  range  p(Torr)  =  0(1-10). 

Since  Av,  /Av,  and  R  are  small  and  large,  respectively,  in  chemical  lasers 
h  d 

the  latter  generally  appear  as  a  product  [e.  g.  ,  Eqs.  (34)]  which  can  be 
approximated  by  (Avh/Avd)  R  =  Ofp(Torr)].  As  previously  noted,  inhomoge¬ 
neous  broadening  effects  become  negligible  in  the  limit  R  — *•  ».  Effects  of 
finite  R  can  be  deduced  from  the  numerical  results  discussed  subsequently. 

Continuous -wave  chemical  laser  performance  has  been  evaluated  for 
laminar  mixing,  a  single  longitudinal  mode  (j^  =1),  0.  1  5  -  0.  5,  0.  01  S  Av^/ 

Av  ,  5  0.  1,  and  10  5  R  <  100.  The  limit  R  -*  00  has  also  been  considered.  The 
d 

results  are  given  in  Table  1  and  Figs.  5  and  6.  Results  for  line-center  opera¬ 
tion,  vij  -»  Vq,  are  given  in  Table  1  and  Figs.  5a,  5b,  6a,  and  6b.  The  approxi¬ 
mate  theory,  i.e.  [l/(l  +  R)]"  «  1,  is  in  good  agreement  with  the  exact  theory 
[Eqs.  (24a)  and  (41)]  for  [l/(1  +  R)]^  <  0(1),  particularly  with  regard  to  net 
output  power  P  .  The  variation  of  the  differential  number  density  2  AN  with 
£  in  the  lasing  region  <  Q  <  is  shown  in  Figs.  5a  and  6a  for  Av^/  Av^  =  0.1 
and  0.01,  respectively.  The  quantity  C^^AN  continues  to  increase,  after 

lasing  is  initiated.  The  increase  is  more  pronounced  for  Av^/  Av^  =  0.01. 

1/2 

This  result  is  in  contrast  with  the  homogeneous  case  (R  -*  ®),where  AN 
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Longitudinal  Mode  j.  =  1  and  ( 


NUMBER  DENSITY 


DISTANCE  C 

Fig.  5a.  Chemical  Laser  Performance  for  Laminar  Mixing, 
Single  Longitudinal  Mode  (jf  =  1),  and  Av^/Av^  =  0.  1 
Equation  (45)  is  used  in  approximate  solution,  i.  e.  , 
[1/(1  +  R)]^  <3C  1.  Equations  (24b)  and  (41)  are  used 
in  exact  solution.  ^  ^AN  versus  Q  for  vj  =  vq. 

5i  £  C  S  Ce- 


DISTANCE  £ 

5b.  Chemical  Laser  Performance  for  Laminar  Mixing, 
Single  Longitudinal  Mode  (j£  =1),  and  Av^/Av^  =  0. 
Equation  (45)  is  used  in  approximate  solution,  i.  e. 
[1/(1  +  R)]^  <?C  1.  Equations  (24b)  and  (41)  are  used 
in  exact  solution.  I  versus  £  for  =  vq. 


POWER  P, 


DISTANCE  C 


Fig.  6a.  Chemical  Laser  Performance  for  Laminar 
Mixing,  Single  Longitudinal  Mode  (j£  =1), 
and  =  0.  01.  Equation  (45)  is  used 

in  approximate  solution,  i.e.,  [1/(1  +  R)]2<Kl 
Equations  (24b)  and  (41)  are  used  in  exact 
solution.  versus  £  for  v>j  =  vg, 

Ci£  CsCe- 


1 


Fig.  6b.  Chemical  Laser  Performance  for  Laminar 
Mixing,  Single  Longitudinal  Mode  (jf  =1), 
and  Av^/Avd  =  0.  01.  Equation  (45)  is  used 
in  approximate  solution,  i.e.,  [ I/(  1  +  R)]^  <$C  1 . 
Equations  (24b)  and  (41)  are  used  in  exact 
solution.  I  versus  £  for  vj  =  vq. 
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FiKil 


FREQUENCY  — ! — si 
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Fig.  6c.  Chemical  Laser  Performance  for  Laminar 

Mixing,  Single  Longitudinal  Mode  (jf  =1),  and 
Av^/Avd  =  0-  01.  Equation  (45)  is  used  in 
approximate  solution,  i.e.,  [1/(1  +  R)]^  1. 

Equations  (24b)  and  (41)  are  used  in  exact 
solution.  Pe  versus  (v  <  =  vq)/Avj. 


remains  constant  in  the  lasing  region.  Similarly,  the  local  lasing  intensity 
I  tends  to  first  increase  with  Q  and  then  decrease  to  zero  (Figs.  5b  and  6b). 

In  the  corresponding  homogeneous  medium,  I  decreases  monatonically  with 
Q  downstream  of  the  region  where  lasing  is  initiated.  The  approximate  theory 
tends  to  underestimate  the  upstream  values  of  I  and  to  overestimate  the  down¬ 
stream  values.  Reasonably  accurate  estimates  of  net  outpower  result 
even  for  those  cases  in  which  the  approximate  theory  for  the  streamwise 
variation  of  I  departs  significantly  from  the  exact  theory.  The  effect  on  out¬ 
put  power  of  operation  off  of  line  center  is  indicated  in  Figs.  5c  and  6c.  For 
R  *  00 ,  the  output  power  has  a  local  minimum  at  because  I+(v^)  and 

I~(Vq)  compete  for  th  '  same  molecules.  The  output  power  first  increases 
with  an  increase  in  since  I+(Vj)  and  I  (v^)  tend  to  interact  with  different 

molecules.  The  power  is  at  maximum  at  v,  -  v.  =  OfAv.  ).  Thereafter,  P 

1  U  n  e 

decreases  with  increase  in  because  of  the  decrease  in  p(v).  The  resultant 
dip  is  similar  to  that  deduced  by  Lamb.  ^  For  a  given  value  of  G  ,  the 
departure  from  the  curves  R  =  00  in  Figs.  5c  and  6c  represents  the  effect  of 
inhomogeneous  broadening  on  output  power.  The  departure  is  small  when 
(Av^/Avj)R  >  10.  Thus,  the  departure  should  be  small  when  cw  chemical 
lasers  operate  in  the  pressure  regime  p(Torr)  >0(10).  Inhomogeneous 
broadening  effects  become  important  for  single -mode  lasers  operating  in  the 
regime  p(Torr)  SO(l).  In  this  pressure  regime,  however,  the  use  of  multiple 
longitudinal  modes  [i.e.,  Av^/Av^  «  1  and  Av^/Av^  =0(1)]  can  result  in 
efficient  power  extraction. 
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The  Lamb  dip  has  been  observed,  experimentally,  in  cw  chemical 
lasers.  The  variation  of  single  mode  output  power  with  lasing  frequency 
is  reproduced  in  Fig.  7  and  is  similar  to  the  curves  in  Figs.  5c  and  6c. 

The  streamwise  variation  of  index  of  refraction  can  also  be  deter¬ 
mined.  Consider  cases  wherein  there  is  no  mode  competition,  i.  e.  , 

Av  »  Av,  .  The  index  of  refraction  at  each  lasing  transition  v.  can  be 
s  h  J 

expressed  [e.  g.  ,  as  Eqs.  (28),  (29c),  and  Appendix  B], 


( 2tt / X)  [ T)  ( v. )  -  1  ]  0.  D(X) 


1/2 

TT  p 


[1  +  0(  Y)] 


where  X,  Y  and  <p.  are  defined  in  Appendix  B  and  g  is  the  threshold  gain 
J  c 

[Eq.  (15)].  Equation  (49)  is  similar  in  form  to  the  result  deduced  from  Lamb's 


theory*^'  wherein  <p.  g  /p.  represents  line  center-zero  power  gain.  The 


local  value  of  index  of  refraction  in  the  present  model  is  seen  to  depend  on 


0.  which  is  a  measure  of  the  local  degree  of  saturation  and  is  evaluated 
from  the  solution  of  Eqs.  (24b)  and  (28).  For  cases  where  (Av^/Av.)  R  »  1, 
Eq.  (48c)  indicates  that  ^  =  1  +  0  j  [  (Avh/Avd)R]~  *{ .  For  the  latter  cases, 
"hole  burning"  effects  are  negligible  since  the  local  index  of  refraction 


equals  the  value  for  an  inhomogeneously  broadened  medium  with  line  center 


gain  equal  to  g  /p y  That  is,  the  index  of  refraction  is  proportional  to  the 
threshold  gain  and  is  unaffected  by  the  zero  power  gain. 
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IV.  CONCLUDING  REMARKS 


It  has  been  assumed  that  the  streamwise  length  of  the  Fabry-Perot 

resonator  equals  the  streamwise  length  of  the  positive  gain  region  in  order 

to  extract  as  much  power  as  possible  from  the  lasing  medium.  In  these 

cases  C  is  of  order  one  (e.  g.  .  Table  1)  and  the  resonator  length  x  is  of 
e  e 

the  order  of  the  characteristic  collisional  deactivation  distance  Xcd  =  u/kcd* 

Thus,  the  particle  transit  time  xg/u  is  of  the  order  of  the  characteristic 

collisional  deactivation  time  k  The  result  that  inhomogeneous  broadening 

cd 

effects  are  small  for  (Av,  /Av  ,)  RiO  (10)  is  therefore  physically  realistic 

n  u 

since  Av,  /Av,  is  a  measure  of  the  fraction  of  the  excited  particles  which 
h  d 

are  resonant  with  the  laser  radiation  field  and  R  is  a  measure  of  the  number 
of  times  the  resonant  and  nonresonant  excited  particles  are  interchanged  by 
cross  relaxation  during  transit  through  the  resonator.  Parameter  values  in 
the  range  (Avh/Avd)  R  >  1  imply  that  all  excited  particles  are  resonant  with 
the  laser  radiation  field  at  some  time  during  transit  through  the  resonator. 

Inhomogeneous  broadening  effects  become  more  severe  in  low  pres¬ 
sure  cw  chemical  lasers  when  the  streamwise  length  of  the  resonator  is 
smaller  than  the  streamwise  length  of  the  positive  gain  region.  In  the  latter 
case  the  number  of  particle  collisions  per  transit  through  the  resonator  is 

of  order  (Ax  /x  ,)  R,  where  Ax  is  the  resonator  length,  and  inhomogeneous 
c  cd  c 

broadening  effects  are  expected  to  be  negligible  when 

(Axc/xcd)  (Av^/Avj)  R>  0(10).  (50) 
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The  latter  reduces  to  the  previous  expression  (Av^/Av^)  R  >  O  (10)  when 

the  resonator  length  is  of  the  order  of  the  streamwise  length  of  the  positive 

gain  region  (i.  e.  ,  when  Ax^  is  of  order  x^).  Eq.  (50)  may  be  viewed  as 

a  generalized  criteria  which  accounts  for  arbitrary  resonator  length.  The 

influence  of  resonator  length  on  inhomogeneous  effects  in  cw  chemical  lasers 

has  been  observed  experimentally.  In  particular,  longitudinal  mode  pulling* 

was  found  to  be  more  severe  when  a  stable  resonator  (Ax  =  0.  1  cm)  was 

c 

1 4 

applied  to  a  cw  chemical  laser  supersonic  gain  medium  (x  ^  =  1  cm)  than 
when  an  unstable  resonator  (Ax^  =  1  cm)  was  applied. 


Chodzko ,  R.  A.  and  Wang,  C.  P.  ,  The  Aerospace  Corporation,  private 
communication. 
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APPENDIX  A 


CHARACTERISTIC  VALUES  OF  PARAMETERS 


Characteristic  values  of  line  shape  parameters,  rate  parameters  and 
stimulated  emission  coefficients  are  given  herein  with  emphasis  on  HF  lasers. 
A.  LINE  SHAPE 
1.  DOPPLER  BROADENING 

The  characteristic  width  of  a  Doppler-broadened  line  shape  is  (FWHM) 


Avd 


=  7.  163  X  10 


[  sec 


-L 


( A  -  1 ) 


where  T  is  in  K,  and  M  is  the  molecular  weight  of  the  lasing  molecule  in 
grams  per  mole.  For  an  HF  laser  (M  =  20,  \  =  2.  77  X  10”^  m),  Av^  =  300  X 
106(T/300)1/2. 

2.  PRESSURE  BROADENING 

The  characteristic  width  of  a  pressure -broadened  line  shape  can  be 
3  1 6 

expressed  ’  (FWHM)  as 


Av^  =  5.996  x  1010£  p.  (atm)  y .  (atm  *  cm  S  (sec  * )  (A-2) 

i 


where  p^  is  the  partial  pressure  of  species  i  and  y^  is  the  broadening  coefficient 
for  that  species.  The  coefficient  y..  depends  on  temperature  level,  as  well  as 
rotational  and  vibrational  energy  levels  J  and  v,  respectively.  Values  of  y^ 
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Table  A-l.  P-Branch  (v  +  1,  J  -  1  -*  v,  J)  Pressure-Broadening  Coefficients 


1.  COLLISIONAL  DEACTIVATION  RATE 

From  Ref.  12,  the  HF  collisional  deactivation  process  is  characterized 
by  the  rate  for  the  process 


r] 

HF(v)  +  HF-^  HF(v  -  1)  +  HF  (A-5a) 

d[Hgv)]  =  -kcd[HF(v)]  (A  -5b) 


,  1Z 

where 


cd 


PHF  “  82.  O^T 


/  3  X  1( 
\  T 


.14 


+  3. 5  X  104T2'  26J 


(A  -6 ) 


The  choice  v  =  2  provides  a  typical  value  for  k 
2.  HARD-SPHERE  COLLISION  RATE 

Let  k  ^  denote  the  number  of  collisions  per  unit  time  that  a  particle  of 
species  A  undergoes  with  the  particles  of  species  B.  The  quantity  (k  *  is 
the  characteristic  time  between  collisions.  For  a  hard  sphere  model, 


gk  _  8.40  X  1 


o25  ( 


°A+° 


B 


/  1 


M 


B 


1/2 


,  -1  .  -1, 
(sec  atm 


(A  -7) 


where  k  ^  is  termed  the  "gas  kinetic"  collision  rate.  Here,  0=1  and  2  are 
for  A  =  B  and  A  *  B,  respectively,  and  a  is  the  particle  diameter  in 
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centimeters.  The  gas  kinetic  rate  for  collisions  between  an  HF  molecule 

-  8 

and  other  HF  molecules  is  (note  that  o  =  2.  7  x  10  cm,  and  0=1) 

Hr  ~ 


Isii  _ 


=  4.47  x  10 


HF 


w 


1/2 


,  "I  t  -1, 

(sec  atm  ) 


(A  -8) 


Collisions  of  HF  with  other  species  are  not  considered  herein.  It  can  be 
assumed  that  the  number  density  distribution  function  n(v)  is  randomized  after 
each  collision  between  HF  and  a  particle  of  similar  or  greater  mass.  Hence, 
assume  =  k^.  Corresponding  values  of  R,  for  an  HF  laser,  are  given 
in  Table  A-2.  It  is  seen  that  10  <  R  <  200  for  100  <  T  <  1000.  Despite  the 
supersonic  expansion  in  cw  chemical  lasers,  the  temperature  in  the  lasing 
region  tends  to  be  of  order  500  K  because  of  reaction  heating  and  nozzle  wall 
boundary  layer  effects.  Thus,  typically,  R  =  0(100). 

3.  STIMULATED  EMISSION  COEFFICIENTS 

The  gain  for  a  P-branch  transition  (v  +  1,  J  -  1  ■*  v,  J)  can  be  written 
for  the  case  of  rotational  equilibrium. 


_  =  a  J  n  ,  - 
r,  J  V,  J\  v  +  1 


-2JT  /T 
e  n 


(A -9) 


An  approximate  expression  for  the  line 

.  12  , 

is  (with  «  Av^  assumed) 


center  value  of  a  T  for  an  HF  laser 

v,  J 


v,  J 


4. 26  X  10 

ZTTT 


11 


(1  +  V  -  0.  Olv  )(1  +  0.  063J)J 
exp [J(J  -  1)TR/T] 


( A  -  1 0) 


ss  Relaxation  Rate  R  =  kcr/kcd  for  HF  Based  on 
•  (A -7)  and  (A -5)  with  v  =  2 


where  =  30.  16  K. 


As  previously  noted. 


Comparison  with  the  present  development  indicates 


o  p  Av, 
o  o  h 


(A-  1 1 ) 


a  Av  is  independent  of  pressure. 
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APPENDIX  B 


MATHEMATICAL  EXPRESSIONS 


The  mathematical  expressions  given  here  have  been  used  in  the  present 
study.  Integrals  involving  the  Lorentzian  line  shape  have  the  form 


00 

I :n  J- 


dv  _  TV 

Av,  -  2 

h 


(B- la) 


(B-  lb) 


00 

f  n*  -  v)'/'(v"-  v)^l  = 


i  +  t(v"-  v')2/(Avh)2] 


(B-lc) 


/CO 

T 


p(v)  <£(\>.  -  v)  dv/Av, 

- __J _ _ _ 

+  [21/(1  +  R)]  ^(Vj  -  v) 


1  (  -X"  2 


-  Y  [1  -  2XD(X)]  +  O  (yV  (B  -  Id) 


4  /*°°  v.  -  v  p(v)^(v.  -  v)  dv/A\^ 

Avh  1  +[2I./(1  +  R)]  <2?(v.  -  v) 


=  ,=  D(X)  -  2XYe'X  +  O  (Y2) 


(B-le) 


where 


X  =  Z  >J  In  2  (vj  -  v^/Av^ 
Y  =  \/  In  2  (Av^/Av^)  <p j 


and  D(X)  is  the  Dawson  integral  which  has  the  following  properties: 


D(z) 


,  /-z  2 

=  e  J  e  dz0 
o 

r  2  2  4  4  ] 

=  n1  -  t  z  +  7t  z  +  •••] 


1  dD[(A;)1/Z]  _ - ‘  -D[(Al)‘/2] 

A  d£ 


2(AC) 


/.£ 

A  J  D[(AC0)1/2]d£0  =  (AC)1/2  -  D[(A4)1/2] 


(B-2a) 

( B  -  2b ) 

( B  -  2c ) 

(B-3a) 

(B  -  3b) 
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NOMENCLATURE 


coefficients,  Eq.  (34) 
speed  of  light  in  vacuum 
Dawson  integral.  Appendix  B 
full  width  at  half  maximum 
normalized  gain,  Eq.  (20) 

gain  per  unit  length  at  frequency  ,  Eq.  (7a) 

local  lasing  intensity  in  ±y  direction 

_  + 

net  local  lasing  intensity,  £  I  (v.) 

j  J 

saturation  intensity,  [(1  +  R)/2]  ek 
normalized  local  intensity,  Eq.  (20) 
longitudinal  mode  number 
final  value  of  j 
characteristic  rate,  sec 

characteristic  collisional  deactivation  rate,  sec 

characteristic  cross  relaxation  rate,  sec 

separation  between  mirrors 

Lorentzian  distribution,  Eqs.  (2b)  and  (20d) 

one  half  and  one  for  laminar  and  turbulent  diffusion, 
respectively 

normalized  number  density  per  unit  volume 

N2-Nl 

particles  (moles)  per  unit  volume  in  interval  v 
to  v  +  dv 


00 


p  (C) 
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e 
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p(vj)-  Pj 
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x,  y 


yf  (*) 
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TKv) 
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characteristic  number  density,  [ F] 
nl  +  n2 

net  particles  (moles)  per  unit  volume  in  upper  and 
lower  lasing  levels 

net  output  power  per  unit  nozzle  height  up  to  station  Q 
net  output  power  per  unit  nozzle  height 
normalized  net  output  power,  Eq.  (29b) 

Maxwellian  distribution  function,  Eq.  (5) 
normalized  Maxwellian  distribution  function,  Eq.  (5) 
static  pressure  in  lasing  region 


k  /k  , 
cr  cd 

mirror  reflectivity 

static  temperature,  °K 

flow  velocity  in  +x  direction 

random  particle  velocity  in  +y  direction 

streamwise  and  transverse  directions,  respectively 

characteristic  diffusion  distance,  Eq.  (17) 

flame  sheet  shape,  Eq.  (17) 

energy  per  mole  of  photons 

normalized  streamwise  distance,  k  ,  x/u 

cd 

normalized  diffusion  distance,  k  ,  x~/u 

cd  D 

index  of  refraction,  Eq.  (7) 
wavelength 


frequency,  sec 
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SUBSCRIPTS 
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00 
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lasing  frequency,  Eqs.  (11)  and  (12) 
longitudinal  mode  separation,  c/2L 
Doppler  width  (FWHM) 
homogeneous  width  (FWHM) 

separation  between  adjacent  modes,  Eq.  (14b) 
cross  section  for  stimulated  emission,  Eq.  (1 

value  at  Vq 

lower  lasing  level 

upper  lasing  level 

value  upstream  of  flame  sheet 

optical  cavity  value 

value  at  end  of  lasing  region 

value  at  start  of  lasing  or  species 

value  corresponding  to  v- 
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LABORATORY  OPERATIONS 

The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver¬ 
satility  and  flexibility  have  been  developed  to  a  high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel¬ 
opments  is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans¬ 
fer,  reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo¬ 
sensitive  materials  and  sensors,  high  precision  laser  ranging,  and  the  appli¬ 
cation  of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optic s ;  communication  sciences,  applied  electronics,  semi¬ 
conducting,  superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor¬ 
rosion  and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia¬ 
tion  from  the  atmosphere,  density  am.  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia¬ 
tions  in  space  on  space  systems. 
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